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ABSTRACT 


Bistatic radar has some properties that are distinetly different from monostatie 
radar. Reeently bistatie radar has reeeived attention for its potential to deteet stealth 
targets due to enhanced target forward scatter. Furthermore, the feasibility of hitchhiker 
radar has been demonstrated, whieh allows passive radar reeeivers to deteet and track 
targets. This thesis developed a software simulation paekage in MATLAB that provides 
a eonvenient tool to examine the bistatie radar design parameters and prediet system 
performanee. The software model is suitable for instruetional purposes due to its user- 
friendly graphieal user interfaee. Several bistatie radar applieations were used to 
illustrate the software features, and their results were analyzed and diseussed. 
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EXECUTIVE SUMMARY 


Bistatic radar uses different antennas separated at a substantial distanee, for 
transmission and reeeption. This fundamental differenee in the bistatie radar has 
generated several properties that are distinctly different from the monostatie radar. These 
unique properties have offered bistatic radar certain advantages in several eommereial 
and military applieations. 

Reeently bistatic radar has received attention for its potential to deteet stealth 
targets due to enhaneed target forward seatter. Furthermore, the feasibility of hitehhiker 
radar has been demonstrated, whieh allows passive radar reeeivers to deteet and traek 
targets, using transmitter of opportunity. Thus, it is likely that bistatic radar will play a 
much bigger role in near-future military warfare. 

Design, development, and testing of radar systems have traditionally been 
expensive and time-eonsuming tasks. Radar designers are often faeed with many design 
iterations to optimize the system performance. Thus, using eomputer simulation and 
modeling will provide the designer the essential tools to examine design parameters and 
to evaluate the system performance before the aetual implementation in hardware. 

There are many software applieations available for monostatie radar simulation 
but none exits in the publie domain for bistatie radar. The objective of this thesis is to 
develop a general software model to analyze the bistatie radar system performance as a 
funetion of the system design parameters. 

The thesis foeused on the design and development of a Graphieal User Interfaee 
(GUIs) in MATLAB to integrate the bistatie radar parameters and to examine several 
eritieal performanee tradeoffs. The end produet is a program named Bidarsa, whieh 
provides a eonvenient instruetional tool for engineering eommunities to understand, 
analyze, and predict bistatie radar performance. 

The thesis presents an overview of bistatie radar definitions and the parameters 
involved in the range equation. The radar deteetion eontour, the ealeulation to determine 
target loeation, eonsiderations for the bistatie footprint, the doppler relationship for a 
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moving target, and the eleetronic countermeasures (ECM) issues are discussed with their 
supporting mathematical equations. These equations were used to formulate the software 
model. 

The Bidarsa program consists of eight major modules developed and integrated to 
carry out the required performance prediction and system analyses. A set of radar 
parameters was compiled to demonstrate the software functions. The thesis further 
illustrates the software functions with two bistatic radar applications. The first 
application used bistatic radar to monitor shipping channel, and the second application 
investigated the feasibility of using direct television broadcast satellites and ground-based 
bistatic receivers to detect targets. 



I. INTRODUCTION 


Bistatic radar has come a long way from its initial days well before the outbreak 
of World War II. Many seientists were trying new radar eonfigurations where the 
reeeiver was located separately from the transmitter. The main motivation at that time 
was to gain sufficient radio frequency isolation between the two electronics sets, 
partieularly with the numerous ongoing eontinuous wave (CW) experiments [1]. This 
fundamental difference in the bistatie radar over monostatic radar has offered many 
opportunities in both commercial and military applications. Over the years and through 
three resurgences, many bistatic radar systems have been proposed, developed, and tested 
in many eountries including the United States, United Kingdom, Franee, Soviet Union, 
Germany and Japan. 

This unique difference has offered bistatic radar certain advantages for particular 
tasks [2]; 

• The receiver is eompletely passive, and hence undeteetable, by the eleetronic support 
measures, and it is safe from attaek by anti-radiation missiles or deliberate directional 
interference and jamming. 

• The reeeiver ean be located in a favorable area where transmitters are not allowed, 
sueh as flammable-liquids stores, gas terminals, ete. 

• No transmit-receive switch or duplexer is required. These devices are lossy, 
expensive, and heavy (important for air-home radar applications). 

• With certain configurations, less transmitter power is required compared to the 
monostatic radar. 

• Higher pulse repetition frequencies (PRFs) can be used because a bistatic system does 
not suffer the same range blindness as the equivalent monostatie system. 

• If the target angle ean be measured at both sites, as well as the bistatic range, data ean 
be cheeked for self-consisteney to remove false alarms. 
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Taking a closer examination of the military applieations, bistatie radar prineiples 
have been deployed in systems sueh as semiaetive homing missiles, forward-seatter 
fences, multistatie radar, and hitchhikers [3]. Mueh research is being done, and 
eontinues, in partieular with regard to the last two applieations mentioned. 

The multistatie radar uses several distributed reeeivers assoeiated with a single or 
multiple transmitters. It is suitable for wide area traeking and it also inereases the 
probability of deteetion. The bistatic hitchhiker uses a transmitter of opportunity, either 
friendly (cooperative) or hostile (non-eooperative), to deteet and loeate targets near the 
transmitter or reeeiver site. There is no doubt that these bistatie radar prineiples are 
playing a bigger role in today’s military for the deteetion of low signature targets and 
improving the eleetronie eountermeasures and eounter-eountermeasures. 

Design, development, and testing of radar systems have traditionally been 
expensive and time eonsuming tasks. Radar designers are often faced with many design 
iterations to optimize the system performanee. Using eomputer simulation and modeling 
provides the designer with the essential tools to examine design parameters and evaluate 
the system performanee before the aetual implementation in hardware. 

There are many software applications available for monostatie radar simulation 
but none exits in the publie domain for bistatie radar. The objeetive of this thesis was to 
develop a general software model to analyze the bistatie radar system performanee as a 
funetion of the system design parameters. 

The thesis foeused on the design and development of a Graphieal User Interfaee 
(GUIs) in MATLAB to integrate the bistatie radar parameters, and examine several 
eritieal performanee tradeoffs. The end produet is a program named Bidarsa, whieh 
provides a eonvenient instructional tool for engineering eommunities to understand, 
analyze, and prediet bistatie radar performanee. 

Chapter II provides the theoretieal baekground required in the formulation of the 
software model. It presents an overview of bistatie radar definitions and the parameters 
involved in the range equation. The radar deteetion eontour, the ealculation to determine 
target loeation, eonsiderations for the bistatie footprint, the doppler relationship for a 
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moving target, and the eleetronic countermeasures (ECM) issues are discussed with their 
supporting mathematical equations. 

Chapter III presents the design approach for the Bidarsa program. A system 
block diagram detailing the interfaces and data flow sequences is used to illustrate the 
functions of the eight major modules and how they are integrated. 

Chapter IV uses a set of radar system parameters compiled to demonstrate the 
Bidarsa program functions and discuss its results. 

Chapter V uses the Bidarsa software to predict the system performance in 
different practical scenarios suggested in References [2] and [4]. The result is analyzed 
and compared with the monostatic radar configuration. 

Chapter VI summarizes the research and suggests further work that could be 
performed to enhance the code. 

Finally, the appendix lists the software codes written for the Bidarsa program. 
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II. BISTATIC RADAR THEORY 


A, DEFINITIONS 

The IEEE Standard 686-1997 [5] has defined bistatic radar as a radar using 
antennas for transmission and reception at sufficiently different locations that the angles 
or ranges from those locations to the target are significantly different, figure 2.1 is used 
to illustrate this definition, showing a transmitter Tx and receiver Rx being situated at two 
locations with a baseline separation distance L. The bistatic angle, P, is one of the 
important parameters which characterizes bistatic radar and affects system performance. 
These relationships will be discussed in the subsequent sections. 



figure 2.1 Bistatic radar geometry. 


There is no strict definition specifying the distance between the transmitter and 
receiver sites. Several discussions [1, 3, 6] suggest that the baseline, L, could range from 
tens of meters to hundreds of kilometers, depending on the required functions of the 
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particular bistatic radar system. The distanee Rj is the range from the transmitter-to- 
target, while Rr is the range from the target-to-receiver. 

B, RANGE EQUATION 

The bistatie range equation is derived in a similar manner as the monostatie radar 
range equation. The bistatie maximum range produet is given by [3] 

P,G,G,l^c,F}Fl 1 ’ p 

or 

)„„=*• ( 2 . 2 ) 

where 

Rt = transmitter-to-target range (m), 

Rr = target-to-reeeiver range (m), 

Pt = transmitted power output (W), 

Gt = transmitting antenna power gain, 

Gr = reeeiving antenna power gain, 

X = transmitting frequeney wavelength (m), 

Ob = target bistatie radar eross seetion (RCS) (m ), 

Ft = pattern propagation factor for the transmitter-to-target path, 

Fr = pattern propagation faetor for the target-to-reeeiver path, 

k = Boltzmann’s eonstant, 1.3807 x 10'^^ J/K, 

Ts = reeeiver noise temperature (°K), 

Bn = reeeiver noise bandwidth (Hz), 

{S/N)min = minimum signal-to-noise power ratio required for deteetion. 
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Lt = transmitting system losses (> 1), 

Lr = receiving system losses (> 1), 

K = bistatic maximum range product. 

In the bistatic radar range equation, it is observed that the maximum range {RfRid^ 
has replaced the range R!* in the monostatic radar equation, due to the different locations 
of the transmitter and receiver. 

The pattern-propagation factors for the transmitter-to-target and target-to-receiver. 
Ft and Fr, also need to be taken care of separately in the bistatic case. This is because 
the transmitting path and receiving path can be very different in nature as compared to 
the monostatic implementation. 


C. OVALS OF CASSINI 


The oval of Cassini is the locus of a points, the product of whose distance from 
two fixed points is a constant. In the bistatic radar case, it is interesting to analyze the 
system performance by plotting the ovals of Cassini as a function of signal-to-noise 
power ratio, S/N, with the ranges, Rt and Rr, at the vertices of the triangle. 

Equation (2.1) is modified and written as 


{RtRr) 


2 

max 


K 




(2.3) 


where K is the bistatic radar constant and given by 

(44kT,B,LrL. 


(2.4) 


The constant K is related to the bistatic maximum range product, K, by 

K = (2.5) 

Figure 2.2 shows Rj and Rr being converted into the polar coordinate system 
{r,d). From geometry 
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Rj^ = + 1 } ! A^ + rL cos 0 


( 2 . 6 ) 


and 


+ 1 } !A^-rLcosG. 


(2.7) 


Combining Equations (2.6) and (2.7) yields 

[R^R^f =(r^ ! A)^ -r^L^ cos'^ 0. 


( 2 . 8 ) 



Figure 2.2 


Bistatic radar geometry converted to the polar coordinate system. (After 

Ref [3].) 


The constant SIN ratio can subsequently be plotted as ovals of Cassini using the 
equation 


(5/A) =-^-. 

(r^+L^IA) -r^L^cos^G 


(2.9) 


Reference [3] has provided an example of the oval of Cassini plot, as shown in Figure 
2.3, with K arbitrarily set to 30C'^ and constant SIN ratios from 10 dB to 30 dB. 
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For analytical purposes, erg, Ft and Fr are assumed to be invariant with r and 6 so 
as to generate a constant-SNR oval-of-Cassini plot. A series of plots shows that bistatic 
radar operation can be divided into three distinct regions; 

1. L> 2-sIlc . Two separate ellipses result enclosing the transmitter and 

receiver. It is considered receiver centered whenR^. transmitter 

centered wheni?^ » Rj.. 

2. L< 2-s[k . A single continuous ellipse results. 

3. L = 2yfK . A lemniscate with a cusp at origin results. 


SVN • to <1B 



Figure 2.3 Ovals of Cassini, with constant SNR plots for K=30L‘^. (From Ref. [3].) 


D. TARGET LOCATION 

A bistatic receiver usually measures the angles of arrival of the target echo and 
estimates the range. The receiver-to-target range, Rr, eannot be measured direetly, but 
can be calculated by solving the bistatic triangle. When the transmitter location is made 
known to the receiver, the baseline range, L, is known. The total signal propagation time 
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from transmitter to target to receiver is also required to be established. This is needed for 
the receiver to convert the time measurement into the range sum estimation, Rt + Rr. 

From Figures 2.2 and 2.3, it can be seen that the isorange contour has to be an 
ellipse, rather than the tangent approximation. It is defined by 

Rj,+Rj^=2ci (2.10) 

where 2a is the major axis of the ellipse, the eccentricity of the ellipse, e, is 

e = L/2a (2.11) 

e = L/{R^ + Rj^). (2.12) 

Using the law of cosines, Rr and Rt can be calculated by measuring the L, (Rt + Rr), and 
either the angle 6t or Or. 

For the case when Or is measured 

Rl=I^+L^- 2R^L cos (90° - 0j ) 

^ (Rt+rJ-l- 

^ 2[Rj +Rj^-LsmOj) 

1 

Rj^ = [rI + 1} -2Rj.L sin . 

Combining Equations (2.12) and (2.14) yields 

^ 2e(l-esin^j,) 

and combining Equations (2.12) and (2.15) yields 

L{^e^ + \-2R sin^^) 

^ 2e[\-esmOj) 

Similarly, for the case when Or is measured 

Rl=Rl+L^- 2R^L cos (90° + 0 ^) 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 
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2[Rj, + Rj^+LsinOi^) 


(2.19) 


R 


Rj = + 1 } + 2Rj^L sin Oj^Y . 

Combining Equations (2.12) and (2.19) yields 

* 2e(l + esin^^) 

Finally, combining Equations (2.12) and (2.20) yields 

+ l + 2i?sin^^) 

^ 2e(l + esin0^) 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


Using the law of sines, it is also possible to establish the parameter ratios to 
calculate the bistatic angle, P, with the measured Rr and Rt and the calculated Or and Or. 
Their relationship is given by 


Thus 


Also, 


Rj^ Rj, L 

sin(9O°-0j,) sin (90° + ^^) sm/3 


(2.23) 



CO?,0j 

(2.24) 

Rj 

cosOg^ 


cosOj 

(2.25) 

L 

sm/3 

sin“^ 

^ L cos 0j ^ 

(2.26) 

\ J 

Rj. 

cosGg^ 

(2.27) 

L ~ 

sin P 
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/5 = sm ^ 


(2.28) 


^ L cos 0^ ^ 

y Rt ) 


E. COVERAGE AREA 

There are two eonsiderations for establishing the eoverage area of the bistatie 
radar: (1) the deteetion-eonstraint whieh is determined by the maximum range oval of 
Cassini, and (2) the line-of-sight (LOS) eonstraint whieh is determined by the geometry 
between the transmitter, target, and reeeiver. 


1. Detection-Constrained Coverage 

The eoverage is determined by the area within a maximum-range oval of Cassini, 
whieh is established by the bistatie maximum range produet, /r, in a bistatie radar 
system. 

Depending on region of operation, i.e., the relationship between L and 24k: , 
different expressions for the area, Ab, apply. The detailed derivation of eaeh expression 
is provided in Appendix D of [3]. 

a. L< 24k: . A single eontinuous ellipse. 


A,,=7tK 


1 - 


^ivr H 




\6k^ 


1 

viy 


1»3 

V2.4y 


r 


\6k^ 


rn 

UJ 


(2.29) 


Ag^ ~ TTK 


1- 

f 1 ^ 





164, 

y ) 


U6384j 



(2.30) 


b. L> 244 . Two separate ellipses enelosing the transmitter and reeeiver. 


^B1 


27CK^ 


1} 


1 + 


-I- 


Af 


2^.2! 


16r" 

IT 


( 



2 


[2^41.3! y 


\ ^ J 


-I- 

J V 

3 


2".3!*2! 


14 


-!-• 


(2.31) 
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(2.32) 


A 


iTtK?' 


B2 




\ 2 k ^ 12r^ lOOr" 

1 + — 


V 


n 


n 
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c. L = I'Jk . A lemniscate with a cusp at the origin 

A^,=2k. (2.33) 

Comparing a bistatic radar with maximum range product, /r, to a monostatic 
radar with an equivalent monostatic range, Vr , it is found that the area of coverage for a 
bistatic radar is smaller than that of the monostatic radar [3]. 


2. Line-of-Sight-Constrained Coverage 

The bistatic line-of-sight (LOS) coverage is affected by multipath, diffraction, 
refraction, and shadowing, including the earth curvature. The coverage area, Ac, is 
established by the common area covered by both the transmitter and receiver coverage 
circles, as shown in Figure 2.4. The 4/3 earth model gives radii of coverage circles as 
follows 

=130(7^ + 7^) (2.34) 





(2.35) 


where 

ht = target altitude (km), 

hr = transmitter antenna altitude (km), 

Hr = receiver antenna altitude (km). 

For L + Tj. > rj^> L-Tj. or L + rj^> Tj, > L-Tj^, the intersection of the two coverage 
circles determines the common coverage area 

= I [<I>R - sin (l>R ) + r^[(l)T- sin (2.36) 
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where 


^J=l cos 


-1 


A 2 1 , 

Vj. +L 

ItjL 


( 1 2 , 

, - _i r„ -Tj. +L 

= 2 cos ^ 


V 


2r^Z 


y 


otherwise, 


if rj. + rj^<L, then = 0; 


(2.37) 

(2.38) 


(2.39) 


if rj.>L + rjj, then ; 


(2.40) 


if , then = TTr^. 


(2.41) 



Figure 2.4 Coverage area geometry. (From Ref [3].) 

F. CUTTER CELL AREA 

The intersection of the range cell and the mainbeam footprint is used to determine 
the bistatic mainlobe clutter cell area, Ac. It is useful to determine if a range cell is 
beamwidth limited or pulsewidth limited. 
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Figure 2.5 illustrates the beamwidth limited case with clutter cell area given by 


[Rj^A9r)[RtA0t) ^2 42 ) 

sm/3 

The terms andi?^A^^ are the cross-range dimension of the transmitting and 

receiving beam at the clutter cell, respectively. The angles AO^ and A9r are the 3-dB 
beamwidths of the transmitting and receiving antennas, respectively. 




Figure 2.5 Beamwidth limited clutter cell area. (From Ref. [7].) 


Figure 2.6 illustrates the pulsewidth limited case with clutter cell area given by 




COS 


(A/ 2 ) 


ctRrAOr 

2cos^ {^/2) 


(2.43) 


The term t is the compressed pulsewidth, whereas AR is the approximation to the 
bistatic range cell given by 
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CT 


(2.44) 




2cos(;5/2) 



TX 


Figure 2.5 Pulsewidth limited clutter cell area. (From Ref. [7].) 


Appendix B of Reference [3] has provided an estimation for the maximum error 
of a range cell, . The maximum error occurs on the perpendicular bisector of the 
baseline, where is a maximum. It is given by 


a{a -a) 
~ b{b'-b) 


(2.45) 


where 


a = semimajor axis of inner ellipse = {Rj + 7?^ ) / 2, 
a = semimajor axis of outer ellipse = (a + cf/ 2), 

b = semiminor axis of inner ellipse = {a^ -l} ! Ay , 

(af-fV4? 


b' = semiminor axis of outer ellipse = 
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G. MAXIMUM UNAMBIGUOUS RANGE AND PRF 

The maximum unambiguous range for a bistatic radar is given by 


[Rj, + Rj^)^ 


c 

PRF 


and the maximum unambiguous PRF is given by 



c 

{Rj.+Rj^)^^ 



c 

+ 2/r(l + coSy5) 


(2.46) 


(2.47) 

(2.48) 


H. DOPPLER RELATIONSHIP 


The doppler shift, fs, is the time rate of change of the total path length of the 
scattered signal, normalized by the wavelength. For bistatic radar, the doppler shift is 
given by 





(2.49) 


/.= 


1 

I 


dRr 

—-+ 

dt 


dRg^ 

dt 


(2.50) 


Figure 2.6 shows a bistatic doppler geometry. The target’s velocity vector has magnitude 
V and aspect angle d. They are being projected to the transmitter and receiver as Vj, Vr 
and dr, Sr respectively. A stationary transmitter and receiver will have the projection of 
the target velocity vector onto the transmitter-to-target LOS and receiver-to-target LOS 
as follows 


^ = Vcos{d-/5l2) (2.51) 


-^ = Vcos{d + pl2). (2.52) 

dt 
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Combining Equations (2.50), (2.51) and (2.52) yields 

f^ = {V! X)[co^{d-l3l2) + co^[d+pl2)] (2.53) 

=(2E//l)cos(^cos(y5/2). (2.54) 


I. ELECTRONIC COUNTERMEASURES 


Bistatic radar can also be subjected to intentional interference initiated by hostile 
transmissions from other sources. Thus it is essential to examine how the performance 
will be affected under such an environment. When a noise jammer is employed, the radar 
receiver will be subjected to the jamming power spectrum density, Jo, given by 




(2.55) 


where 


Pj = jammer transmitted power (W), 
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Gj = j ammer transmit antenna power gain, 

Gr = radar reeeiving antenna power gain, 

Fj = j ammer pattern propagation faetor = Fjfj fj^j , 

Fj = propagation faetor, 

fj = jammer antenna patter faetor, refereneed to the reeeiver site, 

/rj = radar reeeiving antenna patter faetor, refereneed to jammer site, 

Bj = j ammer bandwidth (Hz), 

Rj = j ammer-to-reeeiver range (m), 

Lj = jammer system losses (> 1). 

Assuming the ease of a single jammer, substituting Equation (2.55) into (2.1), the bistatie 
bum-through equation is given by 


{RtRr ) ,, 


PG T B 






F 


LtLr 


4;r{S/N) 


(2.56) 


J. SUMMARY 

The bistatie radar equation is derived in a manner similar to the monostatie radar 
range equation, with the bistatie range (RtFrY replaeing the range R'^ in the monostatie 
radar ease. A eonstant SNR plot has the shape of an oval of Cassini. By resolving the 
bistatie triangle, it is possible to ealeulate the target loeation. The radar eoverage is 
defined by the region where the target is deteetable by the reeeiver and within the LOS of 
both transmitter and reeeiver. The bistatie footprint ean be either beamwidth or 
pulsewidth limited. Both doppler relationships and eleetronie eountermeasures 
eonsiderations were also diseussed. Chapter III uses the theory diseussed in this ehapter 
to establish the approaeh for the bistatie radar simulation, and the integration of the 
various software modules. 
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III. SOFTWARE DESIGN APPROACH 


A, OVERVIEW 

Bidarsa is a bistatic radar simulation developed using the graphical user interfaces 
(GUIs) available in MATLAB. The program is designed to provide easy user interfaees 
for radar designers, engineers, and students to examine the various design parameters for 
a bistatie radar and analyze its performance. 

The bistatic radar model bloek diagram is shown in Figure 3.1. Eight main 
modules were developed to earry out different bistatie radar system analyses as diseussed 
in Chapter II. Figure 3.1 also indieates the data flow sequenees. It illustrates how certain 
modules will need to be exeeuted prior to the others, as their results are used in 
subsequent analyses by the other modules. 

The Bidarsa.m module allows a user to seleet one of the three group funetions: (1) 
specify the system parameters, (2) carry out a basic system performance analysis, or (3) 
to earry out an extended system performanee analysis. 

The Defpara.m module provides a convenient way for a user to speeify the 
bistatie radar’s transmitter and receiver design parameters. It also allows target 
information and the environmental information to be specified. The parameters can be 
updated at any point in time to allow a new scenario to be examined, and its effect on the 
overall system performance evaluated. 

The RngCal.m module eomputes the bistatie maximum range produet, bistatie 
radar constant, and the maximum range-unambiguous pulse repetition frequency (PRF). 
This module allows the user to carry out system analysis on four different design 
parameters: (1) transmitter power, (2) system frequency, (3) target RCS, and (4) signal- 
to-noise ratio. The user is also able to plot the deteetion contour and examine the oval of 
Cassini with four different eonstant SNR plots. 
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Figure 3.1 Bistatic radar system analysis model block diagram. 
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The TgtLoc.m module uses the range-sum value and angle 6 to eompute the 
loeation of the target. It provides the transmitter-to-target and target-to-reeeiver ranges 
and ealeulates the bistatie angle. 

The CovAre.m module evaluates the eoverage areas for both deteotion-eonstrained 
and LOS-eonstrained oases. The module also provides an analysis on the transmitter- 
target-reoeiver height relationship requirement, and the effeots of different baseline 
ranges. 

The FotPrt.m module uses the transmitter-target-reoeiver ranges, antenna 
beamwidth and transmitter pulsewidth, to evaluate the bistatie footprint and olutter area. 
The module examines the results and determines if the deteotion is beamwidth or 
pulsewidth limited. The range oell size and maximum range errors are also oomputed in 
this module. 

The DopRel.m module examines the doppler shift effeots of a moving target. It 
also provides an analysis of the target-velooity aspeot angle versus bistatie angle, and 
how they affeot the doppler frequenoy shift. 

The EWEff.m module provides a means to evaluate the bistatie radar performanoe 
under different jammer environments. The module provides an analysis on the bum- 
through range for different target RCS and jammer power soenarios. 

A Help menu is provided in all modules to assist the user to navigate and 
understand how to use eaoh module. 

Print and Close funotions are also designed into all modules to allow the user to 
make a hard eopy and exit the seleeted menu whenever desired. 

The detailed MATLAB programs for the eight modules are listed in the Appendix 
of this thesis. 
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B. 


BIDARSA MAIN MENU 


To run the Bidarsa GUI, the user has to start MATLAB and set the path of 
direetory to where the Bidarsa eode is stored. Typing Bidarsa at the eommand window 
will exeeute the GUI as shown in Figure 3.2. Among the three groupings, the user should 
define the system parameters first before running the other analysis funetions. An 
explanation of the Bidarsa GUI may be obtained by elieking on the Help button. 



Figure 3.2 main menu GUI. 
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c. 


DESIGN PARAMETERS MENU 


This module is designed to collect the system parameters, target information and 
environmental data. These inputs are generally required by all the other modules. It 
functions as a central database for the program, and converts each parameter to a suitable 
unit for subsequent computations. A pop up window is designed to inform the user that 
new parameters are being loaded for any follow-on analysis. Figure 3.3 shows the GUI 
layout designed for the module. A typical value for each parameter was chosen as the 
default value for demonstration and initial calculation purposes. 



Figure 3.3 Define parameter GUI. 
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D. RANGE CALCULATION MENU 


This module is designed to earry out the most fundamental analysis for a bistatic 
radar system. Several other modules use its result to carry out further analysis. Figure 
3.4 shows the GUI layout designed for the module. With the input data from the 
Defpara.m module, it computes the following: 

1. Bistatic maximum range product, {RjRn , as detailed in Equation (2.1). 

2. Bistatic radar constant, K, as detailed in Equation (2.4). 

3. Maximum range-unambiguous PRE, as detailed in Equation (2.48). 

4. System Analysis for 

a. Transmitting power selection 

b. System frequency selection 

c. Target RCS variation 

d. Receiver SNR variation 

5. Detection contour plot for ovals of Cassini, as detailed in Equation (2.9). 



Eigure 3.4 Range calculation GUI. 
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E. 


TARGET LOCATION MENU 


This module is designed to estimate the location of the target using the range-sum 
value and measured angle 0. Figure 3.5 shows the GUI layout designed for the module. 
Only one of the angles, 6t or Or, is required for the reference, as discussed in Chapter II 
Section D. It computes: 


1. Transmitter-to-target range, as detailed in Equations (2.16) and (2.22). 

2. Target-to-receiver range, as detailed in Equations (2.17) and (2.21). 

3. Bistatic beta angle, as detailed in Equations (2.26) and (2.28). 



Eigure3.5 Target location GUI. 
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F. COVERAGE AREA MENU 

This module is designed to evaluate the region or area on the bistatie plane where 
the target is detectable by the receiver and within the LOS of both transmitter and 
receiver. The user can specify a different altitude for the target, transmitter antenna, and 
receiver antenna for this analysis. Figure 3.6 shows the GUI layout designed for the 
module. It computes; 

1. The region of operation, based on the relationship of L and 2-Jk, as 
discussed in Chapter II Section E. 

2. Coverage area for the detection-constrained case, as detailed in Equations 
(2.30), (2.32), and (2.33). 

3. Coverage area for the EOS-constrained case, as detailed in Equations 
(2.36), (2.39), (2.40), and (2.41). 

4. System analysis for 

a. Transmitter antenna altitude requirements. 

b. Receiver antenna altitude requirements. 

c. Target altitude requirements. 

d. Baseline distance variations. 
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Figure 3.6 Coverage area GUI. 
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G. 


BISTATIC FOOTPRINT AND CLUTTER AREA MENU 


This module is designed to evaluate the bistatic footprint and clutter area of the 
bistatic radar system. The user can specify the transmitter-target-receiver ranges, antenna 
beamwidth and transmitter pulse information for different analysis. Figure 3.7 shows the 
GUI layout designed for the module. It computes: 

1. Whether the system is a beamwidth or pulsewidth limited case. 

2. The bistatic footprint area, as detailed in Equations (2.42) and (2.43). 

3. The range cell size, as detailed in Equation (2.44). 

4. Maximum error in the range cell, as detailed in Equation (2.45). 



Figure 3.7 Clutter area GUI. 
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H, DOPPLER RELATIONSHIP MENU 


This module is designed to evaluate the doppler shift effects of a moving target. 
The user can provide different target velocities and aspect angles for the analysis. Figure 
3.8 shows the GUI layout designed for the module. It computes: 

1. Bistatic doppler shift, as detailed in Equation (2.54). 

2. System analysis for 

a. Target velocity aspect angle variation. 

b. Bistatic angle variation. 

c. Bistatic doppler shift variation. 



Figure 3.8 Doppler relation GUI. 
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I. EW EFFECTS MENU 

This module is designed to evaluate the bistatic radar’s performance when 
subjected to a jamming transmission. The user can define the jammer specifications for 
simulating different operational scenarios. Figure 3.9 shows the GUI layout designed for 
the module. It computes: 

1. Bistatic maximum range product with jammer, as detailed in Equation 

(2.56). 

2. System analysis for 

a. Jammer power variation. 

b. Target RCS variation. 

c. Burn-through range analysis. 



Figure 3.9 EW Effects GUI. 
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J. 


SUMMARY 


The MATLAB GUI is used to develop the Bidarsa simulation program, as it 
offers easy user interfaees and excellent data visualization. Eight modules were 
developed and integrated as a complete suite. The Bidarsa program provides a user the 
options to specify the bistatic radar system parameters, the target information and 
environmental data. Several critical system performance measures are calculated, and 
data presentation tailored to the study of design tradeoffs. The bistatic doppler 
relationship and electronic countermeasures considerations are also included to provide 
extended analyzes. Chapter III provides sample analyzes on a set of radar parameters. 
Plots generated from the different modules will also be discussed. 
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IV. DESCRIPTION OF SOFTWARE OPERATION 


A, TEST DATA 

To evaluate the Bidarsa program, a set of system parameters were eompiled from 
various published sourees [1,3] and used for a detailed system analysis. The test values 
are as shown in Table 1. 


Table 


Transmitting power 

50dBW 

Transmitter antenna gain 

30 dB 

Frequeney 

10 GHz 

Reeeiver min SNR 

0 dBm 

Reeeiver antenna gain 

30 dB 

Reeeiver noise temperature 

300 K 

Reeeiver noise bandwidth 

60 MHz 

Baseline range 

100 km 

Target RCS 

10 dBm^ 

Bistatie beta angle 

1 degree 

Transmitter system loss 

1 dB 

Receiver system loss 

1 dB 

Tx-Tgt pattern propagation factor 

-3.5 dB 

Tgt-Rx pattern propagation factor 

-3.5 dB 


Test parameters used to evaluate Bidarsa program. 


35 



The Bidarsa program is run and the above values are entered and updated in the Defpara 
menu. The results of all module caleulations are displayed in the MATLAB Command 
Window. 

B, BISTATIC RANGE PRODUCT ANALYSIS 

Selecting the Range Calculation menu, the program computed the following 
results for the three range related parameters, as shown in Figure 4.1. Selecting one of 
the four system analysis buttons on the menu will generate the respective plot as shown 
in Figure 4.2. When the plot detection contour button is selected, the program generates 
a plot for the oval of Cassini for constant SNR at 93 dB, 96 dB, 98.65 dB and 102 dB, as 
shown in Figure 4.3. 


Bidarsa - Range Calculation: 

Bistatic Maximum Range Product RtRr is 6.7706e+006 km^. 

Bistatic Radar constant K is 4.5841e+016. 

Maximum Range-Unambiguous PRF is 1.8299e+000 kHz. 

Figure 4.1 Range calculation menu output. 

From Equations (2.48) and (2.46), it is interesting to note that in a bistatic radar 
system, increasing the bistatic angle, ji, will increase the maximum range-unambiguous 
PRF, which will in turn, reduce the bistatic radar maximum unambiguous range. 

Figure 4.3 resembles Figure 2.3 in Reference [3]. With the system parameters 
used in this analysis, the lemniscate {L = 2'Jk) occurs when the SNR is 98.65 dB; two 
separate ellipses enclosing the transmitter and receiver {L> 2Vr ) are observed at a SNR 

of 102 dB; two single continuous ellipses (C < 2Vr) are observed at SNRs of 93 dB and 
96 dB. 
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Figure 4.2 System analysis; transmit power (top left), frequency (top right), target 
RCS (bottom left) and SNR (bottom right). 


Reference [3] also suggests a few possible military applications for each of the 
operating regions. 

1. For Z > 2^/^, with a small oval around the receiver, i.e., the receiver- 
centered region, the system can be used for 

a. Air-to-ground operations using a stand-off transmitter and 
penetrating aircraft with silent receiver. 

b. A semiactive homing missile. 
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2. For L > 2yjlc , with a small oval around the transmitter, i.e., the 
transmitter-centered region, the system can be used for 

a. Monitoring of activity near a non-cooperative transmitter. 

3. For L < 2\[k , with a single continuous ellipse, the system can be used for 

a. Medium-range air defense. 

b. Missile tracking from ground-based sites. 



Figure 4.3 Constant SNR oval of Cassini plot. 
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C. TARGET LOCATION ANALYSIS 

In the Target Location menu, test values for the range-sum and angle Or as shown in 
Table 2, are used for the evaluation. The program eomputed the results for the three 
target loeation related parameters, as shown in Figure 4.4. 


Range-sum -I- ) 

350 km 

Angle Or 

30 degrees 


Table 2. Test values used for target loeation analysis. 


Bidarsa - Target Loeation: 


Theta R Seleeted 


Distanee of Target to Transmitter 

is 2.0938e+002 km. 

Distanee of Target to Reeeiver 

is 1.4063e+002 km. 

The Beta Angle ealeulated 

is 24.43 deg. 

Figure 4.4 

Target loeation menu output. 


The user ean toggle the angle 6 entry and verify from the results that the isorange 
eontour has to be an ellipse, rather than the tangent approximation. The possibility of 
using either Or or Or for eomputation offers more flexibility in the analysis. 
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D, 


COVERAGE AREA ANALYSIS 


In the Coverage Area menu, the test values for the transmitter, reeeiver, and target 
altitude are shown in Table 3. The program eomputed the following results in 
determining the region of operation and the eoverage areas, as shown in Figure 4.5. A 
system analysis of the transmitter-reeeiver-target height requirement versus baseline 
range is also generated as shown in Figure 4.6. 


Table 3. 


Target height 

30 km 

Transmitter antenna height 

10 km 

Reeeiver antenna height 

10 km 


Test values used for eoverage area ana^ 


ysis. 


Bidarsa - Coverage Area: 


Target is in the Single Cosite Region. 


2 

The Deteetion-eonstrained eoverage area is 2.1271e+007 km . 

Rx and Tx is in LOS. 


The LOS-eonstrained eoverage area 

is 3.7384e+006 km^. 


Figure 4.5 Coverage area menu output. 


In this analysis. Figure 4.6 shows that, for a baseline range T of 100 km, an 
adequate LOS is established for a target height > 700 m, only if the transmitting or 
reeeiving antenna is > 100 m. If the baseline is deereased to 50 km (0.5L), an adequate 
LOS is established for a target height > 700 m, when the antenna is near zero altitude. 
However, if the range is inereased to 1200 km (L2L), to aehieve an adequate LOS for 
target height > 700 m, the antenna will need to be raised to at least 180 m. 


40 





Figure 4.6 Transmitter-target-rec eiver height analysis. 
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E, 


BISTATIC FOOTPRINT ANALYSIS 


In the Bistatic Footprint & Clutter Area menu, the radar parameters as shown in Table 4 
are used for the evaluation. The program determines if the bistatic footprint is limited by 
beamwidth or pulsewidth, as shown in Figure 4.7. The range cell size and maximum 
error were also generated. 


Receiver-target range 

100 km 

Transmitter-to-target range 

50 km 

Receiver antenna beamwidth 

10 degrees 

Transmitter antenna beamwidth 

10 degrees 

Transmitter pulsewidth 

30 psec 


Table 4. Test parameters for bistatic footprint analysis. 



Figure 4.7 Bistatic footprint menu output. 


It can be gathered from the Equations (2.42) to (2.45) that the radar pulsewidth 
plays a bigger role than the beamwidth in determining whether the bistatic footprint is 
limited by beamwidth or pulsewidth. This effect can be verified using this module. 
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It is also verified in this module that variation in the antenna beamwidth will only 
affeet the beamwidth limited bistatie footprint area; it has a small effeet on the range eell 
size and does not affeet the maximum error ealeulated. 

F. DOPPLER RELATION ANALYSIS 

In the Doppler Relation menu, the test values shown in Table 5 are used for the 
evaluation. The program eomputes the bistatie doppler shift as shown in Figure 4.8. A 
system analysis for doppler frequeney shift versus different target aspeet angle and 
bistatie angle is generated as shown in Figure 4.9. 


Target veloeity 

300 m/s 

Target veloeity aspeet angle 

80 degrees 


Table 5. Test parameters used for doppler relationship analysis. 


Bidarsa - Doppler Relation: 

Bistatie Doppler Shift is 3.4597e+001 kHz. 

Figure 4.8 Doppler Relation menu output. 

An analysis of Figure 4.9 illustrates the following observations: 

1. For all bistatie angles, P, if the target veloeity aspeet angle, 6, is between 
-90° and +90°, there will be a positive bistatie target doppler. 

2. For all bistatie angles, P, if the target aspeet angle is normal to the bistatie 
biseetor, 6 = ±90°, there will be zero bistatie target doppler. 

3. For all bistatie angles, P < 120°, if the target aspeet angle is eollinear with 
the bistatie biseetor, a maximum bistatie target doppler will oeeur when 
^ = 0° and a minimum will oeeur when d = 180°. 
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Figure 4.9 Doppler frequency shift analysis with different target aspect angle and 

bistatic angle. 
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G. ELECTRONIC COUTNERMEASURES ANALYSIS 

In the EWEffects menu, the jammer specifleations shown in Table 6 are used for 
the evaluation. The program computes the bistatic maximum range product with 
jamming as shown in Figure 4.10. A system analysis for the burn-through ranges with 
different jammer transmit powers and target RCS is generated in Figure 4.11. 


Jammer transmitted power 

30dBW 

Jammer antenna gain 

5 dB 

Jammer propagation factor 

1 dB 

Jammer bandwidth 

1000 MHz 

Jammer-receiver range 

100 km 

Jammer system loss 

3 dB 


Table 6. Test parameters used for electronic countermeasures analysis 


Bidarsa - EW Effects: 

Bistatic Max Range Product with Jamming (RtRrf is 3.6418e+008 km^. 

Eigure 4.10 EW Effects menu output. 
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Figure 4.11 Burn-through range analysis with different target RCS. 


H, SUMMARY 

A set of radar parameters was used to demonstrate the simulation and analysis 
capability provided by the Bidarsa program. The simulation results from the different 
modules were discussed along with their respective plots. The Bidarsa program offers 
flexibility for the user to study how the individual parameters will affect the overall 
system performance. Chapter V presents more examples to illustrate how the Bidarsa 
program can be used to study specific bistatic configurations and scenarios. 
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V. SAMPLE SYSTEMS ANALYSIS 


A, OVERVIEW 

Chapter IV has shown the various components of the software simulation and 
analysis offered by the Bidarsa program. In this chapter, the use of the Bidarsa program 
is further illustrated with two specific bistatic radar applications. 

B, MONITORING OF SHIPPING CHANNEL 

Chapter 13 of Reference [3] has provided a scenario whereby an aircraft uses a 
bistatic radar system to monitor the shipping across a channel. A 1-GHz, 1-W 
transmitting buoy is dropped in a shipping channel such that ships must pass within 1 km 
of it. The aircraft is flying directly overhead listening for the echoes with a receiving 
antenna gain of 20 dB. The buoy has a unity antenna gain and the aircraft receiver 
requires a detection signal strength of -130 dBW. The bistatic RCS of a ship is 30 
dBsm, and all the losses and multipath effects are neglected. 

The above parameters were entered into the Bidarsa program for the analysis. 
The results are tabulated as shown in Table 7. 


Bistatic maximum range product, RtRy 

6.7359 km^ 

Bistatic radar constant, K 

4.5372e-006 

Distance of target-to-transmitter 

1 km 

Distance of target-to-receiver 

6.7359 km 

Beta angle 

82.17 degrees 

Table 7. Results for shipping-channel 

. example. 


The above bistatic radar result was used to compare its performance to an aircraft 
using monostatic radar in the same scenario. The transmitting antenna is assumed to 
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have a unity gain for consistency. The monostatic radar will give a maximum detection 
range of only 2.595 km, due to the factor of in the monostatic radar range equation. 
The bistatic radar system in this application has a range advantage over the monostatic 
radar. 

This comparison has implied that by using the bistatic radar configuration, the 
aircraft could fly at a much higher altitude to monitor the shipping channel without being 
detected. On the other hand, if the monostatic radar uses a higher gain antenna, it will 
make up for the monostatic factor and require less transmitting power for the same 
range performance. However, the monostatic radar configuration, being an active 
system, might possibly give away the position of the aircraft. 

The transmitting power, system frequency, target RCS, and SNRmin versus 
maximum range product are shown in Figure 5.1. 



Figure 5.1 System analysis for shipping-channel example: transmit power (top left), 
frequency (top right), target RCS (bottom left) and SNR (bottom right). 
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Depending on the requirements of the aircraft operation profile, Figure 5.1 
illustrates the tradeoffs in the system design parameters, so as to achieve a higher 
probability of detection and aircraft survivability. 

C. NFS BISTATIC RADAR EXPERIMENT 

An experiment was conducted at the Naval Postgraduate School (NPS) in 1997 to 
investigate the feasibility of using direct television broadcast satellites and ground-based 
bistatic receivers to detect targets [4]. With the Flughes Direct Broadcast Satellite (DBS- 
1) used to illuminate the target, receive signals from the direct path antenna and bistatic 
receive antenna were compared and processed to detect the presence of targets. The 
system configuration is shown in Figure 5.2. 

The system parameters used for the experiment are as shown in Table 8. The 
parameters were entered into the Bidarsa program for the analysis. The resulting bistatic 
maximum range product is 74543 km , and the estimated RR is approximately 2 km. 
This is in agreement with the result obtained in the experiment conducted. 



Satellite 


Target 


Direct 


Receiver and Processor 


Figure 5.2 NPS bistatic radar experiment configuration. 
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Effective radiated power of DBS-1, PtGt 

53 dBW 

Transmitting frequency,/ 

12.7 GHz 

Effective aperture of receiver antenna, Ae 

0.16 m^ 

Processor gain of radar, Gp 

61 dB 

Receiver bandwidth. Bn 

45 MHz 

Minimum SNR required for detection, (S/N)min 

13 dB 

Distance between transmitter and receiver, L 

37000 km 

Target bistatic RCS, ob 

30 dBsm 

Two way Eosses, E/r 

6dB 


Table 8. System parameters for the NPS bistatie radar experiment. 


The experimental result confirmed the feasibility of the concept, and suggested 
several improvement factors [4] that could employed to increase the detection range. 
They are: 

1. Using of low-noise receiver to improve the SNR. 

2. Increase the receiver antenna area and gain. 

3. Increase the effective radiated power on the transmitting satellite. 

4. Using multiple receivers to increase the processing gain. 

It is important to take note that in the bistatic-radar-satellite-hitchhiker system, each 6-dB 
increase in SNR doubles the receiver-to-target range, Rr. This is much better than the 
monostatic system, whereby it takes a 12-dB increase in SNR to double the detection 
range. 
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D, 


SUMMARY 


In this chapter, two different applications were presented and the Bidarsa program 
was used to carry out the system performance prediction. The results were analyzed and 
eompared with the monostatic radar system. Recommendations were also diseussed to 
improve the system performance. The final chapter will present a overall summary of 
this thesis and suggest possible follow up work. 
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VI. SUMMARY AND CONCLUSIONS 


A, SUMMARY 

Bistatic radar offers many advantages over monostatic radar in certain 
applications, since the reeeiving antenna is at a distant site with reference to the 
transmitting antenna. This thesis achieved the objective of developing a software model 
for bistatic radar systems. The software allows a user to examine the various system 
design parameters and evaluate the system performance tradeoffs, under different 
conditions and scenarios, in a simulated environment. 

MATLAB GUI is chosen as the software to develop the Bidarsa simulation 
program because it is an effective interactive tool for both numerical computation and 
data visualization. It also allows a very user-friendly interface to be designed and 
operated. 

The Bidarsa program is easy to learn and navigate with Help menus designed to 
guide the user in every module. It is suitable to be used by any engineering student or 
professionals who are interested in the bistatic radar field. 

With the eurrent operation trends toward a lower signature target and the third 
resurgence for bistatic radar using Passive Coherent Location [8], the bistatic radar offers 
a much greater potential now than at any time before. The Bidarsa program will help 
engineering students to understand the fundamental design principles of bistatic radar, 
and provide a convenient tool for the professional in this field to analyze the system 
performance. The software is available free of eharge and ean be downloaded from 
www.nps.navv.mil/ienn . 

B, RECOMMENDATIONS 

While this thesis focused on the effort to develop a software model to evaluate the 
essential performanee of a bistatie radar system, there are other areas that eould be 
examined to provide a more advanced analysis to predict the bistatie radar performanee. 
They include: 
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1. Examine the effects of bistatic RCS for surface clutter, which consists of 
both ground and sea echoes. 

2. Clutter tuning concepts that cater to a moving transmitting and receiver. 

3. Pulse chasing concepts to reduce the complexity and cost of multibeam 
bistatic receivers. 

4. Passive Coherent Location (PCL) considerations [8], such as the 
CELLDAR concept shown in Eigure 6.1. CELLDAR is a passive system based on the 
bistatic radar principle. It works by using existing transmissions made by digital cell 
phone base stations. When a target object enters the detection region, these cell phone 
transmissions are reflected by it and detected at the CELLAR antennas. 


1. Aircraft enters 


CElIDAP 



ihi 


Cell Phone 
Radar 


Acoustic 

Detection 

Sensor 


Cell 

Phone 

Base Station 


Central 

Control System 


Eigure 6.1 CELLDAR concept of operation proposed by Roke Manor. (Erom Ref. 

[9].) 
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APPENDIX 


A, MATLAB CODE FOR BIDARSA.M 

This appendix details the Bidarsa.m eode developed to allow a user to seleet one 
of the three group funetions: (1) speeify the system parameters, (2) earry out a basie 
system performanee analysis, or (3) earry out an extended system performanee analysis. 


funetion varargout = Bidarsa(varargin) 

% BIDARSA M-file for Bidarsa.fig 

% BIDARSA, by itself, ereates a new BIDARSA or raises the existing 
% singleton*. 

% H = BIDARSA returns the handle to a new BIDARSA or the handle to 
% the existing singleton*. 

% BIDARSA('CALLBACK',hObjeot,eventData,handles,...) ealls the loeal 
% funetion named CALLBACK in BIDARSA.M with the given input arguments. 
% BIDARSA('Property','Value',.--) creates a new BIDARSA or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before Bidarsa OpeningFunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to Bidarsa OpeningFcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help Bidarsa 
% East Modified by GUIDE v2.5 20-0ct-2003 17:08:59 
% Begin initialization code - DO NOT EDIT 
guiSingleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'gui_Singleton', guiS ingle ton, ... 

'gui_OpeningEcn', @Bidarsa_OpeningEcn, ... 

'gui_OutputEcn', @Bidarsa_OutputEcn, ... 

'gui EayoutEcn', [] , ... 

'gui_Callback', []); 
if nargin & isstr(varargin{ 1}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Lnargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin ):} ); 
end 

% End initialization code - DO NOT EDIT 
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% — Executes just before Bidarsa is made visible. 

function Bidarsa_OpeningFcn(hObject, eventdata, handles, varargin) 

handles, output = hObject; 

guidata(hObject, handles); 

clear all; 

clc; 


% — Outputs from this function are returned to the command line, 
function varargout = Bidarsa_OutputFcn(hObject, eventdata, handles) 
varargoutjl} = handles, output; 


% — Executes on button press in pushbuttonOO. 

function pushbuttonOO_Callback(hObject, eventdata, handles) 

Defpara 

% — Executes on button press in pushbutton 11. 

function pushbuttonl l_Callback(hObject, eventdata, handles) 

RngCal 

% — Executes on button press in pushbuttonl2. 

function pushbutton 12_Callback(hObject, eventdata, handles) 

TgtEoc 

% — Executes on button press in pushbuttonl3. 

function pushbuttonl3_Callback(hObject, eventdata, handles) 

CovAre 

% — Executes on button press in pushbuttonl4. 

function pushbutton 14_Callback(hObject, eventdata, handles) 

FotPrt 

% — Executes on button press in pushbutton21. 

function pushbutton21_Callback(hObject, eventdata, handles) 

DopRel 

% — Executes on button press in pushbutton31. 

function pushbutton31_Callback(hObject, eventdata, handles) 

EWEff 
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% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

close all hidden 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpMain 
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B, 


MATLAB CODE YOKDEFPARA.M 


This appendix details the Defpara.m code developed to provide a convenient way 
for user to specify the bistatic radar’s transmitter and receiver design parameters, the 
target information and the environmental information. 


function varargout = Defpara(varargin) 

% DEFPARA M-file for Defpara.fig 

% DEFPARA, by itself, creates a new DEFPARA or raises the existing 
% singleton*. 

% H = DEFPARA returns the handle to a new DEFPARA or the handle to 
% the existing singleton*. 

% DEFPARA('CAEEBACK',hObject,eventData,handles,...) calls the local 
% function named CAEEBACK in DEFPARA.M with the given input arguments. 
% DEFPARA('Property','Value',...) creates a new DEFPARA or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before Defpara OpeningFunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to Defpara OpeningFcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help Defpara 
% East Modified by GUIDE v2.5 24-Oct-2003 01:25:10 
% Begin initialization code - DO NOT EDIT 
gui_Singleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui_OpeningFcn', @Defpara_OpeningFcn, ... 

'gui_OutputFcn', @Defpara_OutputFcn, ... 

'gui EayoutFcn', [],... 

'gui Callback', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Enargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin):} ); 
end 

% End initialization code - DO NOT EDIT 


% — Executes just before Defpara is made visible. 

function Defpara_OpeningFcn(hObject, eventdata, handles, varargin) 
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handles, output = hObject; 
guidata(hObject, handles); 
ele; 


% — Outputs from this function are returned to the command line, 
function varargout = Defpara_OutputFcn(hObject, eventdata, handles) 
varargoutjl} = handles, output; 


% — Executes during object creation, after setting all properties, 
function Pt_input_CreateFcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Pt_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Gt_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Gt_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Ec[_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Eq_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function SN_input_CreateEcn(hObject, eventdata, handles) 
if ispc 
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set(hObject,'BackgroundColorVwhite'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function SN_input_Callback(hObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Gr_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Gr_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Ts_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Ts_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Bn_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Bn_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Br_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 
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set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Br_input_Callback(hObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Sm_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Sm_input_Callback(bObject, eventdata, bandies) 


function Ba_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultFbcontrolBackgroundColor')); 

end 

function Ba_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Ft_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultFbcontrolBackgroundColor')); 

end 

function Ft_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Fr_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 
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function Fr_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Lt_input_CreateFcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Lt_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Er_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function Er_input_Callback(hObject, eventdata, handles) 


% — Executes on button press in Elpdate. 

function Update_Callback(hObject, eventdata, handles) 

global Pt OO Gt OO Eq_00 SN OO Pt OO Gr OO Ts OO Bn OO Br OO... 

Sm OO Ba OO Ft OO Fr OO Et OO Er OO Eambda OO 
clc; 

Pt_input = eval(get(handles.Pt_input,'String')); 

Gt_input = eval(get(handles.Gt_input,'String')); 

Fq_input = eval(get(handles.Fc[_input,'String')); 

SN_input = eval(get(handles.SN_input,'String')); 

Grinput = eval(get(handles.Gr_input,'String')); 

Ts_input = eval(get(handles.Ts_input,'String')); 

Bninput = eval(get(handles.Bn_input,'String')); 

Br_input = eval(get(handles.Br_input,'String')); 

Sminput = eval(get(handles.Sm_input,'String')); 

Ba_input = eval(get(handles.Ba_input,'String')); 

Ft_input = eval(get(handles.Ft_input,'String')); 

Fr_input = eval(get(handles.Fr_input,'String')); 

Etinput = eval(get(handles.Et_input,'String')); 

Er_input = eval(get(handles.Er_input,'String')); 
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Pt_00 = 10^(Pt_input/10); 

GtOO = 10^(Gt_input/10); 

Fc[_00 = Fq_input*le9; 

SN_00 = (10^((SN_input-30)/10)); 

BaOO = deg2rad(Ba_input); 

GrOO = 10^(Gr_input/10); 

TsOO = Tsinput; 

BnOO = Bn_input*le6; 

BrOO = Brinput; 

Sm_00 = 10^(Sm_input/10); 

Ft_00 = 10^(Ft_input/10); 

FrOO = 10^(Fr_input/10); 

LtOO = 10^(Lt_input/10); 

LrOO = 10^(Lr_input/10); 

Lambda_00 = 3e8/(Fq_00); %System Wavelength 

mbO = msgbox([' New bistatie radar parameters loaded. '],'Result'); 


% — Exeeutes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 

% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpDefpara 


63 



C. MATLAB CODE FOR RNGCAL.M 

This appendix details the Rngcal.m eode developed to compute the bistatic 
maximum range product, bistatic radar constant, maximum range-unambiguous PRF, 
plotting of the detection contour, and carry out system analysis on four design 
parameters. 


function varargout = RngCal(varargin) 

% RNGCAL M-file for RngCal.fig 

% RNGCAL, by itself, creates a new RNGCAL or raises the existing 
% singleton*. 

% H = RNGCAL returns the handle to a new RNGCAL or the handle to 
% the existing singleton*. 

% RNGCAL('CALLBACK',hObject,eventData,handles,...) calls the local 
% function named CALLBACK in RNGCAL.M with the given input arguments. 
% RNGCAL('Property','Value',.--) creates a new RNGCAL or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before RngCal OpeningFunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to RngCal OpeningFcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help RngCal 
% East Modified by GUIDE v2.5 06-Nov-2003 08:56:10 
% Begin initialization code - DO NOT EDIT 
gui_Singleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'gui_Singleton', guiS ingle ton, ... 

'gui_OpeningEcn', @RngCal_OpeningEcn, ... 

'gui_OutputEcn', @RngCal_OutputEcn, ... 

'gui EayoutEcn', [] , ... 

'gui_Callback', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Lnargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin ):} ); 
end 

% End initialization code - DO NOT EDIT 
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% — Executes just before RngCal is made visible. 

function RngCal_OpeningFcn(hObject, eventdata, handles, varargin) 

handles, output = hObject; 

guidata(hObject, handles); 

clc; 

axes(handles. axes 1); 

[xl,map] = imread('z-Range.jpg',jpg'); 

image(xl); 

colonnap(map); 

set(handles.axesl,'VisibleVoff); 
global PowP FreP RCSP SNRP 
PowP = 0; 

FreP = 0; 

RCSP = 0; 

SNRP = 0; 


% — Outputs from this function are returned to the command line, 
function varargout = RngCal_OutputFcn(hObject, eventdata, handles) 
varargoutjl} = handles, output; 


% — Executes on mouse press over figure background, over a disabled or 

% — inactive control, or over an axes background. 

function figurel_WindowButtonDownFcn(hObject, eventdata, handles) 


% — Executes on key press over figure 1 with no controls selected, 
function figurel_KeyPressFcn(hObject, eventdata, handles) 


% — Executes on button press in Calculate. 

function Calculate_Callback(hObject, eventdata, handles) 

global Pt OO Gt OO Fq_00 SN OO Pt OO Gr OO Ts OO Bn OO Br OO... 

Sm OO Ba OO Ft OO Fr OO Et OO Er OO Eambda OO RtRrl 1 BRPl 1... 

PowP FreP RCSP SNRP ConP 
clc; 

k= 1.38065*10^(-23); 

Numll = (Pt_00*Gt_00*Gr_00*Eambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Denll = ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Et_00*Er_00); 

RtRrl 1 = ((Numl 1/Denl 1)^(0.5))/1000; %Bistatic Max Range Product in km 

Denl2 = ((4*pi)^3*k*Ts_00*Bn_00*Et_00*Er_00); 

BRPl 1 = (Numl 1/Denl2); %Bistatic Radar Parameter in km 

PRFull = 3e8/((Br_00^2+(2*(RtRrll*1000)*(l+cos(Ba_00))))^0.5); 

PRFul2 = PRFul 1/1000; 
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fprintf('\nBidarsa - Range CalculationAn'); 

fprintf('Bistatic Maximum Range Produet RtRr is %-8.4e km^2.\n',RtRrl 1); 
fprintf('Bistatic Radar eonstant K is %-8.4e.\n\n',BRPl 1); 
fprintf('Maximum Range-Unambiguous PRF is %-8.4e kHz.\n\n',PRFul2); 


olose(figure(l)); 

%Plotting Transmit Power Analysis 
if (PowP==l) 

Ptll = (Pt_00/10):(Pt_00/10):(Pt_00M0); 

Numl 11 = (Gt_00.*Pt_l 1 *Gr_00*Lambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Benin = ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl 11 = (((Numl 11 ./Deni 11).^(0.5))/1000); %Bistatic Max Range Product in km 

Pt_112 = 10.*logl0(Pt_ll); 

figure(l); 

plot(RtRrlll,Pt_112); 
title('Transmit Power Analysis'); 
xlabel('Max Range Product (km^2)'); 
ylabel('Transmit Power (dBW)'); 
grid on; 
end 

%Plotting Tansmit Frequency Analysis 
if (FreP==l) 

Fq_12 = (Fq_00-5e9):(Fq_00/10):(Fq_00+5e9); 

Lambda_12 = 3e8./Fq_12; 

Numl21 =(Pt_00*Gt_00*Gr_00.*Lambda_12.^2*Sm_00*Ft_00^2*Fr_00^2); 
Denl21 = ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl21 = (((Numl21./Denl21).^(0.5))/1000); %Bistatic Max Range Product in km 

Fq_121 =Fq_12./le9; 

figure(l); 

plot(RtRrl21,Fq_121); 
title('Transmit Frequency Analysis'); 
xlabel('Max Range Product (km^2)'); 
ylabel('Freq (GHz)'); 
grid on; 
end 


%Plotting Target RCS Analysis 
if(RCSP==l) 

Sm_13 = (Sm_00/10):(Sm_00/10):(Sm_00*10); 

Numl31 = (Pt_00*Gt_00*Gr_00*Lambda_00^2.*Sm_13*Ft_00^2*Fr_00^2); 

Denl31 = ((4*pi)^3*k*Ts_00*Bn_00*SN_00*Lt_00*Lr_00); 

RtRrl31 = (((Numl31./Denl31).^(0.5))/1000); %Bistatic Max Range Product in km 
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Sm_121 = 10.*logl0(Sm_13); 
figure(l); 

plot(RtRrl31,Sm_121); 
title('Target RCS Analysis'); 
xlabel('Max Range Product (km^2)'); 
ylabel('Target RCS (dBsm)'); 
grid on; 
end 

%Plotting Receiver SNR Analysis 
if(SNRP==l) 

SN_14 = (SN_00/10):(SN_00/10):(SN_00*10); 

Numl41 = (Pt_00*Gt_00*Gr_00*Lambda_00^2*Sm_00*Ft_00^2*Fr_00^2); 

Denl41 = ((4*pi)^3*k*Ts_00*Bn_00.*SN_14*Lt_00*Lr_00); 

RtRrl41 = (((Numl41./Denl41)A(0.5))/1000); %Bistatic Max Range Product in km 
SN_141 = (10.*logl0(SN_14))+30; %Converting dB to dBm 

figure(l); 

plot(RtRrl41,SN_141); 
title('SNR_{(min)} Analysis'); 
xlabel('Max Range Product (km^2)'); 
ylabel('SNR_{(min)} (dBm)'); 
grid on; 
end 


% — Executes on button press in Pow_plot_input. 

function Pow_plot_input_Callback(hObject, eventdata, handles) 


set(hObject, 'Value', 1); 
set(handles.Fre_plot_input, 'Value', 0); 
set(handles.RCS_plot_input, 'Value', 0); 
set(handles.SNR_plot_input, 'Value', 0); 
global PowP FreP RCSP SNRP 
PowP = 1; 

FreP = 0; 

RCSP = 0; 

SNRP = 0; 


%Turn this button on 
%Turn off all the buttons 
%Turn off all the buttons 
%Turn off all the buttons 


% — Executes on button press in Ere_plot_input. 
function Ere_plot_input_Callback(hObject, eventdata, handles) 
set(hObject, 'Value', 1); %Turn this button on 

set(handles.Pow_plot_input, 'Value', 0); %Turn off all the buttons 

set(handles.RCS_plot_input, 'Value', 0); %Tum off all the buttons 

set(handles.SNR_plot_input, 'Value', 0); %Tum off all the buttons 

global PowP EreP RCSP SNRP 
PowP = 0; %get(handles.Pow_plot_input, 'Value'); 
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FreP = 1; %get(hObject,'Value'); 

RCSP = 0; 

SNRP = 0; 

% — Executes on button press in RCS_plot_input. 
function RCS_plot_input_Callback(bObject, eventdata, bandies) 
set(bObject, 'Value', 1); %Turn this button on 

set(bandles.Pow_plot_input, 'Value', 0); %Turn off all tbe buttons 

set(bandles.Fre_plot_input, 'Value', 0); %Tum off all tbe buttons 

set(handles.SNR_plot_input, 'Value', 0); %Tum off all the buttons 

global PowP FreP RCSP SNRP 
PowP = 0; 

FreP = 0; 

RCSP = 1; 

SNRP = 0; 

% — Executes on button press in SNR_plot_input. 
function SNR_plot_input_Callback(hObject, eventdata, handles) 
set(hObject, 'Value', 1); %Turn this button on 

set(handles.Pow_plot_input, 'Value', 0); %Tum off all the buttons 

set(handles.Fre_plot_input, 'Value', 0); %Tum off all the buttons 

set(handles.RCS_plot_input, 'Value', 0); %Tum off all the buttons 

global PowP FreP RCSP SNRP 
PowP = 0; 

FreP = 0; 

RCSP = 0; 

SNRP= 1; 

% — Executes on button press in Con_plot_input. 

function Con_plot_input_Callback(hObject, eventdata, handles) 

close(figure(2)); 

CalculateCallback; 

Cass 


% plot ovals of cassini for bistatic radar SNR 

function Cass 

global BRPl 1 Br OO 

rad=pi/180; 

i=0; 

F=Br_00; 

K=BRP11; 

delt=.l; 

a0= 10Nogl0(16*K/(FM)); 
al = round(aO); 
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am6 = al-6; 
am3 = al-3; 
ap3 = al+3; 
snr=[am6 am3 aO ap3]; 

for isnr=l:length(snr) 

SNR= 10^(snr(isnr)/10); 
i=i+l; k=0; 

% setup for b~=a 
a=L/2; b4=K/SNR; b=b4^(.25); 


if b>=a 

for thd=0:delt;360 
thr=thd*rad; 

% b>a is one curve 
if b>a 
k=k+l; 

T(k)=thr; 

Rsq=a^2*(cos(2*thr)+sqrt(cos(2*thr)^2-l+(b/a)M)); 

rl(k)=sqrt(Rsq); 

r2(k)=0; 

end 

% b=a is lemniscate 
if b==a, 
k=k+l; 

T(k)=thr; 

Rsq=a^2*2*cos(2*thr); rl(k)=sqrt(Rsq); r2(k)=0; 
end 
end 

figure(2) 

if i==l, hold on; plot(-a,0,'or',a,0,'or'); end 
X=rl.*cos(T); Y=rl.*sin(T); 
plot(X,Y); 
end 
end 

clear rl r2 T X Y 
% b<a is two curves 

% this covers b<a in which case 8 separate curves are generated 
if b<a 

for curv=l;4 

thl=[0 90 180 270]; %lower limit 

th2=[90 180 270 360]; %upper limit 

k=0; 

% decrease angle increment for these small curves 
for thd=th 1 (curv): delt/10; th2(curv) 
thr=thd*rad; 
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Rlsq=a^2*(cos(2*thr)+sqrt(cos(2*thr)^2-l+(b/a)^4)); 

R2sq=a^2*(cos(2*thr)-sqrt(cos(2*thr)^2-l+(b/a)M)); 

Rl=sqrt(Rlsq); R2=sqrt(R2sq); 

if imag(Rl)==0 & imag(R2)==0 %R1~= conj(R2) 

k=k+l; 

T(k)=tbr; 

rl(k)=(Rl); r2(k)=(R2); 
end 
end 

figure(2); X=rl.*eos(T); Y=rl.*sin(T); 

plot(X,Y); 

elear X Y 

if i==l, bold on; end 
X=r2.*eos(T); Y=r2.*sin(T); 
plot(X,Y); 
grid on; 

elear rl r2 T X Y 
end 
end 

title('Detection Contours Analysis - Ovals of Cassini'); 

xlabel('Range (km)'); 

ylabel('Range (km)'); 

axesCposition', [0 0 1 1],'visible','off); 

text(0.31,0.48,'Tx'); 

text(0.70,0.48,'Rx'); 

bold off 

fprintf('Ovals of Cassini illustrate tbe eontour plots for S/N at:\n'); 
fprintf('%-4.2fdB, %-4.2fdB, %-4.2fdB & %-4.2fdB.\n\n',am6, am3, aO, ap3); 


% — Exeeutes on button press in Print. 

funetion Print_Callbaek(bObjeet, eventdata, bandies) 

printdlgC-setup',bandies.figurel) 

% — Exeeutes on button press in Close. 

funetion Close_Callbaok(bObjeot, eventdata, bandies) 

delete(bandles. figure 1) 

elose(figure(l)) 

elose(figure(2)) 

% — Exeeutes on button press in Help. 

funetion Help Callbaek(bObieet, eventdata, bandies) 

yHlpRngCal 
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D, MATLAB CODE FOR TGTLOC.M 

This appendix details the TgtLoc.m eode developed to eompute the loeation of the 

target. 


funetion varargout = TgtLoe(varargin) 

% TGTLOC M-file for TgtLoe.fig 

% TGTLOC, by itself, ereates a new TGTLOC or raises the existing 
% singleton*. 

% H = TGTLOC returns the handle to a new TGTLOC or the handle to 
% the existing singleton*. 

% TGTLOC('CALLBACK',hObjeot,eventData,handles,...) ealls the loeal 
% funetion named CALLBACK in TGTLOC.M with the given input arguments. 
% TGTLOC('Property','Value',...) ereates a new TGTLOC or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before TgtLoe OpeningFunetion gets ealled. An 
% unreeognized property name or invalid value makes property applieation 
% stop. All inputs are passed to TgtLoe OpeningFen via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instanee to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help TgtEoe 
% East Modified by GUIDE v2.5 29-Oet-2003 15:06:14 
% Begin initialization eode - DO NOT EDIT 
guiSingleton = I; 

gui_State = struot('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui OpeningEen', @TgtEoo_OpeningEon, ... 

'gui OutputEen', @TgtEoo_OutputEon, ... 

'gui_EayoutEon', [] , ... 

'gui_Callbaok', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callbaok = str2funo(varargin{l}); 
end 

if nargout 

[varargout)Lnargout}] = gui_mainfon(gui_State, varargin):}); 
else 

gui_mainfen(gui_State, varargin ):} ); 
end 

% End initialization eode - DO NOT EDIT 


% — Exeeutes just before TgtEoe is made visible. 

funetion TgtEoo_OpeningEon(hObjeot, eventdata, handles, varargin) 

handles, output = hObjeet; 
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guidata(hObject, handles); 
global TR TT 
axes(handles. axes 1); 

[xl,map] = imread('z-TgtLoe.jpg','jpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

TR= 1; 

TT = 0; 


% — Outputs from this funetion are returned to the eommand line, 
funetion varargout = TgtLoe_OutputFen(hObjeet, eventdata, handles) 
varargoutjl} = handles, output; 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion BR_input_CreateFon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion BR_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion TR_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion TR_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion TA_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 
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function TA_input_Callback(hObject, eventdata, handles) 


% — Executes on button press in Theta_R_input. 

function Theta_R_input_Callback(hObject, eventdata, handles) 

set(hObject, 'Value', 1); %Turn this button on 

set(handles.Theta_T_input, 'Value', 0); %Turn off all the buttons 

global TR TT 

TR = get(hObject,'Value'); 

TT = get(handles.Theta_T_input, 'Value'); 


% — Executes on button press in Theta T input. 
function Theta_T_input_Callback(hObject, eventdata, handles) 
set(hObject, 'Value', 1); %Turn this button on 

set(handles.Theta_R_input, 'Value', 0); %Turn off all the buttons 

global TR TT 

TR = get(handles.Theta_R_input, 'Value'); 

TT = get(hObject,'Value'); 


% — Executes on button press in Calculate. 

function Calculate_Callback(hObject, eventdata, handles) 

global Br OO TR TT 

clc; 

TR_input = eval(get(handles.TR_input,'String')); 

TA_input = eval(get(handles.TA_input,'String')); 

BR_21 = Br_00; %Base Range in km 

TR_21 = TRinput; 

TA_21 = deg2rad(TA_input); 

Elli =BR_21/TR_21; 

Num21 = BR_21*(1 - E11C2); 

Den21 = 2*Elli*(l + Elli * sin(TA_21)); 

Rr21 = Num21 /Den21; %Tgt Range to Tx in km 

Num22 = BR_21*(E11C2 + 1 + (2*Elli*sin(TA_21))); 

Rt21 = Num22/Den21; %Tgt Range to Rx in km 

Den23 = 2*Elli*(l - Elli * sin(TA_21)); 

Rt22 = Num21/Den23; 

Num24 = BR_21*(E11C2 + 1 - (2*Elli*sin(TA_21))); 

Rr22 = Num24/Den23; 

fprintf('\nBidarsa - Target EocationAn'); 
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ifTR>=TT 
Rt23 = Rt21; 

Rr23 = Rr21; 

Bet_21 = asin(BR_21/Rt23*cos(TA_21)); 
fprintf('Theta_R Selected\n'); 
else 

Rt23 = Rt22; 

Rr23 = Rr22; 

Bet_21 = asin(BR_21/Rr23*cos(TA_21)); 
fprintf('Theta_T Selected\n'); 
end 

fprintf('Distance of Target to Transmitter is %-8.4e km.\n',Rt23); 
fprintf('Distanee of Target to Receiver is %-8.4e km.\n',Rr23); 

Bet_22 = rad2deg(Bet_21); 

fprintf('The Beta Angle calculated is %-0.2f deg.\n\n', Bet_22); 


% — Executes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 

% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpTgtEoc 
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E, MATLAB CODE FOR COVARE.M 

This appendix details the CovAre.m eode developed to evaluate the eoverage 
areas of both detection-constrained and LOS-constrained cases. 


function varargout = CovAre(varargin) 

% COVARE M-file for CovAre.fig 

% COVARE, by itself, creates a new COVARE or raises the existing 
% singleton*. 

% H = COVARE returns the handle to a new COVARE or the handle to 
% the existing singleton*. 

% COVARE('CAEEBACK',hObject,eventData,handles,...) calls the local 
% function named CAEEBACK in COVARE.M with the given input arguments. 
% COVARE('Property','Value',...) creates a new COVARE or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before CovAre OpeningEunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to CovAre OpeningEcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help CovAre 
% East Modified by GUIDE v2.5 23-Oct-2003 11:49:58 
% Begin initialization code - DO NOT EDIT 
guiSingleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui_OpeningEcn', @CovAre_OpeningEcn, ... 

'gui_OutputEcn', @CovAre_OutputEcn, ... 

'gui_EayoutEcn', [] , ... 

'gui_Callback', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Enargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin ):} ); 
end 

% End initialization code - DO NOT EDIT 


% — Executes just before CovAre is made visible. 

function CovAre_OpeningEcn(hObject, eventdata, handles, varargin) 

handles, output = hObject; 
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guidata(hObject, handles); 
axes(handles. axes 1); 

[xl,map] = imread('z-Centerr.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

axes(handles. axes2); 

[x2,map] = imread('z-Centert.jpgVjpg'); 

image(x2); 

eolormap(map); 

set(handles.axes2,'VisibleVoff); 

axes(handles.axes3); 

[x3,map] = imread('z-Centere.jpgVjpg'); 

image(x3); 

eolormap(map); 

set(handles.axes3,'VisibleVoff); 

axes(handles.axes4); 

[x4,map] = imread('z-LOS.jpg',’jpg'); 

image(x4); 

eolormap(map); 

set(handles.axes4,'Visible',’off); 


% — Outputs from this funetion are returned to the eommand line, 
funetion varargout = CovAre_OutputFen(hObjeet, eventdata, handles) 
varargoutjl} = handles, output; 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Htg_input_CreateFon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 
else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion Htg_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Htx_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 
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else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

funetion Htx_input_Callback(hObject, eventdata, bandies) 


% — Exeeutes during object creation, after setting all properties, 
function Hrx_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Hrx_input_Callback(bObject, eventdata, bandies) 


% — Executes on button press in Calculate. 

function Calculate_Callback(bObject, eventdata, bandies) 

global RtRrl 1 Br OO 

clc; 

Htginput = eval(get(bandles.Htg_input,'String')); 
Htx_input = eval(get(bandles.Htx_input,'String')); 
Hrx_input = eval(get(bandles.Hrx_input,'String')); 


Htg_31 = Htginput; 

Htx_31 = Htx_input; 

Hrx_31 = Hrxinput; 
BMRP31 = RtRrl 1; 

BR31 =Br_00; 

TEST = 2*sqrt(BMRP31); 
TEST2 = abs(BR31-TEST); 
TEST3 = TEST2/BR31; 


%Heigbt in km 


%Bistatic Max Range Parameter in km 
%Base Range in km 


fprmtf('\nBidarsa - Coverage Area:\n'); 

if (TEST3<=0.1) %CheckifE~=2sqrt(k) 

Ab31 = pi*BMRP31; %Area in km^2 

fprintfCTarget is in tbe Eemniscate Region.V); 
elseif(BR31<TEST) 

Ab31 = pi*BMRP31*(l - (BR31M/BMRP31^2/64) - (3*BR31^8/BMRP31M/16384)); 
fprintfCTarget is in tbe Single Cosite Region.V); 
else (BR31>TEST) 


77 



Ab31 = (2*pi*BMRP31^2/BR31^2)*(l + (2*BMRP31^2/BR31M) + ... 

(12*BMRP31M/BR31^8) + (100*BMRP31^6/BR31^12)); 
fprintf('Target is in the Two separate Ellipses Region.V); 
end 

fprintf('The Deteetion-eonstrained eoverage area is %-8.4e km^2.\n\n',Ab31); 


Rr31 = 130*(sqrt(Htg_31) + sqrt(Hrx_31)); %Range in km 
Rt31 = 130*(sqrt(Htg_31) + sqrt(Htx_31)); 

Theta_r = 2*aeos((Rr31^2 - Rt31^2 + BR31^2)/(2*Rr31*BR31)); 
ThetaJ = 2*aeos((Rt31^2 - Rr31^2 + BR31^2)/(2*Rt31*BR31)); 

if (Rt31>=(BR31+Rr31)) 

Ae31 = pi*Rr31^2; %Area in km^2 

fprintf('Tx eoverage ineludes all Rx eoverage eirele.W); 
elseif (Rr31 >=(BR31 +Rt31)) 

Ae31 = pi*Rt31^2; 

fprintf('Rx eoverage ineludes all Tx eoverage eirele.V); 
elseif ((Rr31 +Rt31 )>=BR31) %Area in km^2 

Ae31 = 0.5*((Rr31^2*(Theta_r - sin(Theta_r))) + ... 

(Rt31^2*(Theta_t - sin(Theta_t)))); 
fprintf('Rx and Tx is in LOS.V); 
else 

Ae31 = 0 

fprintf('No LOS between Rx and Tx.V) 
end 

fprintf('The LOS-eonstrained coverage area is %-8.4e km^2.\n\n',Ac31); 


close(figure(l)); 

%Plotting Transmitter-Target-Receiver Height Analysis 

BMRP32 = RtRr 11 * 1 e3; %Bistatic Max Range Parameter in m 

BR32 =Br_00*le3; %Base Range in m 

for L= [(0.5*BR32) BR32 (1.2*BR32)] 
n=0; 

for ht=50;5; 1000 %target heigth in m 

n=n+l; 

hrm=( ((BMRP32+L^2/4)^0.5+L/2)/130-sqrt(ht*le3) )^2; 
hr=hrm/1000; 

Ht(n)=ht; 

Hr(n)=hr; 

end 

figure(l) 

plot(Ht,Hr,'-'), hold on 
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title('Target Height Analysis'); 
xlabel('Target Altitude (m)') 
ylabel('Tx or Rx Antenna Altitude (m)') 
grid on 
end 

axesCposition', [0 0 1 1],'visible',’off); 

text(0.30,0.15,'0.5*L'); 

text(0.48,0.20,'L'); 

text(0.68,0.26,'1.2*L'); 

hold off 


% — Executes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 

% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

close(figure(l)) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpCovAre 
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F. MATLAB CODE FOR FOTPRT.M 

This appendix details the ForPrt.m eode developed to evaluate the bistatie 
footprint and elutter area. 


funetion varargout = FotPrt(varargin) 

% FOTPRT M-file for FotPrt.fig 

% FOTPRT, by itself, ereates a new FOTPRT or raises the existing 
% singleton*. 

% H = FOTPRT returns the handle to a new FOTPRT or the handle to 
% the existing singleton*. 

% FOTPRT('CALLBACK',hObjeot,eventData,handles,...) ealls the loeal 
% funetion named CALLBACK in FOTPRT.M with the given input arguments. 
% FOTPRT('Property','Value',...) ereates a new FOTPRT or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before FotPrt OpeningFunetion gets ealled. An 
% unreeognized property name or invalid value makes property applieation 
% stop. All inputs are passed to FotPrt OpeningFen via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instanee to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help EotPrt 
% East Modified by GUIDE v2.5 20-0et-2003 17:10:38 
% Begin initialization eode - DO NOT EDIT 
guiSingleton = I; 

gui_State = struot('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui_OpeningEon', @EotPrt_OpeningEon, ... 

'gui OutputEen', @EotPrt_OutputEon, ... 

'gui_EayoutEon', [] , ... 

'gui_Callbaok', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callbaok = str2funo(varargin{l}); 
end 

if nargout 

[varargout)Lnargout}] = gui_mainfon(gui_State, varargin):}); 
else 

gui_mainfen(gui_State, varargin ):} ); 
end 

% End initialization eode - DO NOT EDIT 


% — Exeeutes just before EotPrt is made visible. 

funetion PotPrt_OpeningPon(hObjeot, eventdata, handles, varargin) 

handles, output = hObjeet; 
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guidata(hObject, handles); 
axes(handles. axes 1); 

[x,map] = imread('z-BWLimit.jpgVjpg'); 

image(x); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 

axes(handles. axes2); 

[y,map] = imread('z-PWLimit.jpgVjpg'); 

image(y); 

eolormap(map); 

set(handles.axes2,'VisibleVoff); 


% — Outputs from this funetion are returned to the eommand line, 
funetion varargout = FotPrt_OutputFen(hObjeet, eventdata, handles) 
varargoutjl} = handles, output; 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Htx_input_CreateFon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion Htx_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Hrx_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion Hrx_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Rr_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 
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set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Rr_input_Callback(hObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Rt_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Rt_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Rx_BW_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Rx_BW_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Tx_BW_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Tx_BW_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function PusWid_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 
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function PusWid_input_Callback(hObject, eventdata, handles) 


% — Executes on button press in Calculate. 

function Calculate_Callback(hObject, eventdata, handles) 

global Ba OO Br OO 

clc; 

Rrxinput = eval(get(handles.Rr_input,'String')); 

Rtxinput = eval(get(handles.Rt_input,'String')); 

RBWinput = eval(get(handles.Rx_BW_input,'String')); 

TBW_input = eval(get(handles.Tx_BW_input,'String')); 

PuW_input = eval(get(handles.PusWid_input,'String')); 

Bet_41 = Ba_00; %Beta Angle in rad 

Rrx_41 = Rrxinput; 

Rtx_41 = Rtxinput; 

RBW_41 = deg2rad(RBW_input); 

TBW_41 = deg2rad(TBW_input); 

PuW_41 = PuW_input* 10^-6; 

Num41 = Rrx_41 *RBW_41 *Rtx_41 *TBW_41; 

Den41 = sin(Bet_41); 

Num42 = 3e8*PuW_41*Rrx_41*RBW_41; 

Den42 = 2*(cos(Bet_41/2))^2; 

DelR = (3e8*PuW_41)/(2*cos(Bet_41/2)); 

RrDr = Rrx_41*RBW_41; 

CeS41 = 3e8*PuW_41/(2*cos(Ba_00/2)); 

Br41 =Br_00*le3; 

Aa41 = (Rrx_41+Rtx_41)*le3/2; 

Ap41 = Aa41 + (3e8*PuW_41/2); 

Bb41 = (Aa41^2 - (Br41^2)/4)^(0.5); 

Bp41 = (Ap41^2 - (Br41^2)/4)^(0.5); 

Emax41 = ((Aa41*(Ap41-Aa41))/(Bb41*(Bp41-Bb41)))-!; 

fprintf('\nBidarsa - Bistatic EootprintAn'); 
if DelR < RrDr 
Ac41 = Num42/Den42; 
fprintf('This is a Pulsewidth Eimited case.\n'); 
else 
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Ac41 = Num41/Den41; 
fprintf('This is a Beamwidth Limited caseAn'); 
end 

fprintf('The Bistatie Footprint is %-8.4e sqr-kmAn\n',Ao41); 

fprintf('Range Cell Size is %-8.4e m.\n',CeS41); 
fprintf('Maximum Error is %-8.4e mAn\n',Emax41); 


% — Executes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 

% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpEotPrt 
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G. MATLAB CODE FOR DOPREL.M 

This appendix details the DopRel.m code developed to examine the doppler shift 
of a moving target. 


function varargout = DopRel(varargin) 

% DOPREL M-fde for DopRel.fig 

% DOPREL, by itself, creates a new DOPREL or raises the existing 
% singleton*. 

% H = DOPREL returns the handle to a new DOPREL or the handle to 
% the existing singleton*. 

% DOPREL('CALLBACK',hObject,eventData,handles,...) calls the local 
% function named CALLBACK in DOPREL.M with the given input arguments. 
% DOPREL('Property','Value',.--) creates a new DOPREL or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before DopRel OpeningLunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to DopRel OpeningLcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDLES 
% Edit the above text to modify the response to help DopRel 
% Last Modified by GUIDE v2.5 02-Nov-2003 11:11:15 
% Begin initialization code - DO NOT EDIT 
gui_Singleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui OpeningEcn', @DopRel_OpeningEcn, ... 

'gui OutputEcn', @DopRel_OutputEcn, ... 

'gui_LayoutEcn', [] , ... 

'gui Callback', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Enargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin):} ); 
end 

% End initialization code - DO NOT EDIT 


% — Executes just before DopRel is made visible. 

function DopRel_OpeningEcn(hObject, eventdata, handles, varargin) 

handles, output = hObject; 
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guidata(hObject, handles); 
axes(handles. axes 1); 

[xl,map] = imread('z-DopRel.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 


% — Outputs from this funetion are returned to the eommand line, 
funetion varargout = DopRel_OutputFen(hObjeet, eventdata, handles) 
varargoutjl} = handles, output; 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Vel_input_CreateFon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultUieontrolBaekgroundColor')); 

end 

funetion Vel_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Del_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion Del_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion TA_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion TA_input_Callbaok(hObjeot, eventdata, handles) 
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% — Executes on button press in Calculate. 

function Calculate_Callback(hObject, eventdata, handles) 

global Ba OO Lambda OO 

clc; 

Vel_input = eval(get(handles.Vel_input,'String')); 

Delinput = eval(get(handles.Del_input,'String')); 

Vel_51 = Velinput; 

Del_51 = deg2rad(Del_input); 

DoS_51 = (2*Vel_51/Lambda_00)*(cos(Del_51)*cos(Ba_00/2)); 
DoS_52 = DoS_51/1000; %in kHz 

fprintf('\nBidarsa - Doppler RelationAn'); 
fprintf('Bistatic Doppler Shift is %-8.4e kHz.\n\n',DoS_52); 


close(figure(l)); 

%Plotting Doppler Shift Analysis 
forBa_51 = 0:(pi/4);pi; 

Del_52 = -pi:(pi/360):pi; 

DoS_53 = (2*Vel_51/Lambda_00).*cos(Del_52).*cos(Ba_51/2); 
figure(l); 

DoS_54 = DoS_53./le3; 

Del_53 = l.*rad2deg(Del_52); 
plot(Del_53,DoS_54); 

title('Bistatic Target Doppler Shift Analysis'); 
xlabel('Target velocity aspect angle, \delta, (deg)'); 
ylabel('Bistatic Doppler Shift, F_{tgt}(kHz)'); 
xhm([-180 180]); 
grid on; 
hold on 
end 

axesCposition', [0 0 1 l],'visible','off); 
text(0.87,0.54,'\beta = 180\circ'); 
text(0.87,0.40,'\beta = 135\circ'); 
text(0.87,0.25,'\beta = 90\circ'); 
text(0.87,0.17,'\beta = 45\circ'); 
text(0.87,0.13,'\beta = 0\circ'); 
hold off 


% — Executes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 
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% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

close(figure(l)) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpDopRel 
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H. MATLAB CODE FOR EWEFF.M 

This appendix details the EWEjf.m code developed to evaluate the bistatic radar 
performance under different jammer environments. 


function varargout = EWEff(varargin) 

% EWEEE M-fde for EWEff fig 

% EWEEE, by itself, creates a new EWEEE or raises the existing 
% singleton*. 

% H = EWEEE returns the handle to a new EWEEE or the handle to 
% the existing singleton*. 

% EWEEE('CAEEBACK',hObject,eventData,handles,...) calls the local 
% function named CAEEBACK in EWEEE.M with the given input arguments. 
% EWEEE('Property','Value',.--) creates a new EWEEE or raises the 
% existing singleton*. Starting from the left, property value pairs are 
% applied to the GUI before EWEff OpeningEunction gets called. An 
% unrecognized property name or invalid value makes property application 
% stop. All inputs are passed to EWEff OpeningEcn via varargin. 

% *See GUI Options on GUIDE'S Tools menu. Choose "GUI allows only one 
% instance to run (singleton)". 

% See also: GUIDE, GUIDATA, GUIHANDEES 
% Edit the above text to modify the response to help EWEff 
% East Modified by GUIDE v2.5 23-Oct-2003 14:10:38 
% Begin initialization code - DO NOT EDIT 
gui_Singleton = I; 

gui_State = struct('gui_Name', mfilename, ... 

'guiSingleton', guiS ingle ton, ... 

'gui_OpeningEcn', @EWEff_OpeningEcn, ... 

'gui_OutputEcn', @EWEff_OutputEcn, ... 

'gui_EayoutEcn', [] , ... 

'gui Callback', []); 
if nargin & isstr(varargin{l}) 

gui_State.gui_Callback = str2func(varargin{l}); 
end 

if nargout 

[varargout)Enargout}] = gui_mainfcn(gui_State, varargin):}); 
else 

gui_mainfcn(gui_State, varargin):} ); 
end 

% End initialization code - DO NOT EDIT 


% — Executes just before EWEff is made visible. 

function EWEff_OpeningEcn(hObject, eventdata, handles, varargin) 

handles, output = hObject; 
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guidata(hObject, handles); 
axes(handles. axes 1); 

[xl,map] = imread('z-EWEff.jpgVjpg'); 

image(xl); 

eolormap(map); 

set(handles.axesl,'VisibleVoff); 


% — Outputs from this funetion are returned to the eommand line, 
funetion varargout = EWEff_OutputEen(hObjeet, eventdata, handles) 
varargoutjl} = handles, output; 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Pj_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColorVwhite'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEheontrolBaekgroundColor')); 

end 

funetion Pj_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Gj_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion Gj_input_Callbaok(hObjeot, eventdata, handles) 


% — Exeeutes during objeet ereation, after setting all properties, 
funetion Ej_input_CreateEon(hObjeot, eventdata, handles) 
if ispe 

set(hObjeet,'BaekgroundColor','white'); 

else 

set(hObjeet,'BaekgroundColor',get(0,'defaultEbeontrolBaekgroundColor')); 

end 

funetion Ej_input_Callbaok(hObjeot, eventdata, handles) 
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% — Executes during object creation, after setting all properties, 
function Bn_input_CreateFcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColorVwhite'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Bn_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Ej_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function Ej_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Rj_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function Rj_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Bj_input_CreateEcn(hObject, eventdata, handles) 
if ispc 

set(hObject,'BackgroundColor','white'); 

else 

set(hObject,'BackgroundColor',get(0,'defaultEhcontrolBackgroundColor')); 

end 

function Bj_input_Callback(hObject, eventdata, handles) 


% — Executes during object creation, after setting all properties, 
function Sm_input_CreateEcn(hObject, eventdata, handles) 
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if ispc 

set(hObject,'BackgroundColorVwhite'); 

else 

set(hObjeet,'BackgroundColor',get(0,'defaultUicontrolBaekgroundColor')); 

end 

function Sm_input_Callback(hObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Lr_input_CreateFcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

function Lr_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Et_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Et_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Er_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 

else 

set(bObject,'BackgroundColor',get(0,'defaultEbcontrolBackgroundColor')); 

end 

function Er_input_Callback(bObject, eventdata, bandies) 


% — Executes during object creation, after setting all properties, 
function Et_input_CreateEcn(bObject, eventdata, bandies) 
if ispc 

set(bObject,'BackgroundColor','wbite'); 
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else 

set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor')); 

end 

funetion Ft_input_Callbaek(hObject, eventdata, bandies) 


% — Exeeutes on button press in Caleulate. 

function Calculate_Callback(bObject, eventdata, bandies) 

global Pt OO Gt OO Fq_00 SN OO Pt OO Gr OO Ts OO Bn OO Sm OO... 

Ft OO Fr OO Ft OO Fr OO Fambda OO RtRrl 1 
clc; 

Pj_input = eval(get(bandles.Pj_input,'String')); 

Gj_input = eval(get(bandles.Gj_input,'String')); 

Fj_input = eval(get(bandles.Fj_input,'String')); 

Bj_input = eval(get(bandles.Bj_input,'String')); 

Rj_input = eval(get(bandles.Rj_input,'String')); 

Fj_input = eval(get(bandles.Fj_input,'String')); 

Pj_61 = 10^(Pj_input/10); 

Gj_61 = 10^(Gj_input/10); 

Fj_61 = 10^(Fj_input/10); 

Bj_61 = Bj_input*10^6; 

Rj_61 = Rj_input*10^3; 

Fj_61 = 10^(Fj_input/10); 

Num61 = Pt_00*Gt_00*(Ft_00)^2*(Fr_00)^2*Fj_61 *Bj_61 *Sm_00*(Rj_61)^2; 

Dem61 = Pj_61*Gj_61*(Fj_61)^2*Ft_00*Fr_00*Bn_00*4*pi*SN_00; 

RtRrJl = (Num61/Dem61)^(0.5); 
fprintf('\nBidarsa - EW Effects :\n'); 

fprintf('Bistatic Max Range Product with Jamming (RtRr)j is %-8.4e km^2.\n\n',RtRrJl); 


close(figure(l)); 

%Plotting EW Effects Analysis 

for Sm_611 = [(Sm OO/lO) Sm OO (Sm_00*10)] 

Pj_611 = (Pj_61/100):(Pj_61/100):(Pj_61*100); 

Num611 = Pt_00*Gt_00*(Ft_00)^2*(Fr_00)^2*Ej_61 *Bj_61 .*Sm_611 *(Rj_61)^2; 
Dem611 = Gj_61.*Pj_61 l*(Fj_61)^2*Et_00*Er_00*Bn_00*4*pi*SN_00; 

RtRrJ61 = ((Num611./Dem611)A(0.5)); 

Pj_612 = 10.*logl0(Pj_611); 
figure(l); 

plot(RtRrJ61,Pj_612); 
title('Burn-tbrougb Range Analysis'); 
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xlabel('Max Range Product with Jamming (km^2)'); 
ylabel('Jammer Transmit Power (dBW)'); 
grid on; 
hold on 
end 

axesCposition', [0 0 1 l],'visibleVoff); 
text(0.15,0.15,'0.1 \sigma_ {B }'); 
text(0.38,0.20,'\sigma_{B}'); 
text(0.60,0.25,'10\sigma_{B}'); 
hold off 


% — Executes on button press in Print. 

function Print_Callback(hObject, eventdata, handles) 

printdlg('-setup',handles.figurel) 

% — Executes on button press in Close. 

function Close_Callback(hObject, eventdata, handles) 

delete(handles. figure 1) 

close(figure(l)) 

% — Executes on button press in Help. 

function Help_Callback(hObject, eventdata, handles) 

yHlpEWEff 
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